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Abstract 

A (2+l)-dimensional quasilinear system is said to be 'integrable' if it can be 
decoupled in infinitely many ways into a pair of compatible n-component one- 
dimensional systems in Riemann invariants. Exact solutions described by these 
reductions, known as nonlinear interactions of planar simple waves, can be viewed 
as natural dispersionless analogs of n-gap solutions. It is demonstrated that the 
requirement of the existence of 'sufficiently many' n-component reductions pro- 
vides the effective classification criterion. As an example of this approach we clas- 
sify integrable (2+l)-dimensional systems of conservation laws possessing a convex 
quadratic entropy. 
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1 Introduction 



In this paper we address the problem of integrability of (2 + 1) -dimensional quasilinear 

systems 

u t + A(u)u x + B(u)uy = (1) 

where t, x, y are independent variables, u is an m-component column vector and A(u), B(u) 
are mxm matrices. Systems of this type describe many physical phenomena. In particu- 
lar, important examples occur in gas dynamics, shallow water theory, combustion theory, 
nonlinear elasticity, magneto-fluid dynamics, etc. |28j . Although many interesting sys- 
tems of the form arise as dispersionless limits of multidimensional soliton equations 
[39J or within the R-matrix approach |3] and, therefore, should be regarded as integrable, 
no intrinsic definition of the integrability for multidimensional quasilinear systems has 
been proposed until recently. In particular, the standard symmetry approach [321 
which proves to be extremely effective in the case of higher order evolution equations and 
systems, does not seem to work in this context. 

The key element of our construction are exact solutions of the system (0) of the form 
u = u(R l ,...,R n ) where the Riemann invariants R 1 ,...,R n solve a pair of commuting 
diagonal systems 

R\ = A*(i?) R x , R l y = n\R) R % x \ (2) 

notice that the number of Riemann invariants is allowed to be arbitrary! Thus, the 
original 2 + 1-dimensional system (JIJ is decoupled into a pair of diagonal 1 + 1-dimensional 
systems. Solutions of this type, known as nonlinear interactions of n planar simple waves, 
were extensively investigated in gas dynamics and magnetohydrodynamics in a series of 
publications [3 0111 HUH EUl HE 120] • Later, they appeared in the context of the dispersionless 
KP hierarchy [THl HE1 H3 HH1 EU EU EE EH! an d the theory of integrable hydrodynamic- 
type chains |3U ES|- We will call a multidimensional system integrable if it possesses 
'sufficiently many' n-component reductions of the form for arbitrary n (the precise 
definition follows). The corresponding nonlinear n-wave interactions can be viewed as 
dispersionless analogs of £ n-gap' solutions. 

We recall, see jHHj, that the requirement of the commutativity of the flows (J2J) is 
equivalent to the following restrictions on their characteristic speeds: 

3 - ^ d^d/d w - (3) 



\' - A' fp - //' ' 

(no summation!). Once these conditions are met, the general solution of (J2J) is given by 
the implicit 'generalized hodograph' formula [3B] 

v%R) =x + X\R) t + fi^R) y, i = l,...,n, (4) 

where v l (R) are characteristic speeds of the general flow commuting with (J2J), that is, the 
general solution of the linear system 

V = d ^ = <v . , . (5) 
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Substituting u(R l , ...,R n ) into (JTJ) and using (j2J), one readily arrives at the equations 

(A + f/B + A*I m ) <9;U = 0, i = l,...,n, (6) 

implying that both X 1 and satisfy the dispersion relation 

det{A + fiB + XI m ) = 0. (7) 

Thus, the construction of nonlinear interactions of n planar simple waves consists of two 
steps: 

(1) Reduce the initial system (JTJ) to a pair of commuting flows (0) by solving the equations 
(JHJ), © for u(R), X l (R), //(-R) as functions of the Riemann invariants R 1 1 ... 1 R n . For 
n > 3 these equations are highly overdetermined and do not possess solutions in general. 
However, once a particular reduction of the form (J2J) is constructed, the second step is 
fairly straightforward: 

(2) Solve the linear system (jSJ) for v l {R) and determine R 1 ,...,R n as functions of t,x,y 
from the implicit hodograph formula (J3J. 

Remark 1. For n = 1 we have u = u(R), where the scalar variable R = R 1 solves a 
pair of first order PDEs 

Rt = A(i?) R x , R y = n(R) Rx 

which, in one-component situation, are automatically commuting. We recall that in the 
scalar case the hodograph formula takes the form 

f(R)=x + \(R)t + ii(R)y (8) 

where f(R) is arbitrary. This formula shows that, in coordinates t,x,y, the surfaces 
R = const are planes. Hence, u = u(R) is constant along a one-parameter family of 
planes. Solutions of this type, known as planar simple waves, exist for all multidimensional 
quasilinear systems and, therefore, cannot detect the integrability. 

Similarly, for n = 2, we have u = u(i? 1 , R 2 ) where R 1 , R 2 satisfy the system (j2J whose 
general solution is given by the generalized hodograph formula 

v 1 (R) = x + X 1 (R)t + fi 1 (R)y, v 2 (R) = x + X 2 (R)t + /j 2 (R)y, R=(R\R 2 ). (9) 

Setting R 1 = const, R 2 = const, one obtains a two-parameter family of lines (or, in 
the geometric language, a line congruence) in the 3-space of independent variables t, x, y. 
Therefore, the solution u = u(i? 1 ,i? 2 ) is constant along the lines of a two-parameter 
family. The requirement of the existence of solutions of this type, known as nonlinear 
interactions of two planar simple waves, is also not very restrictive. For instance, for 
m = 2, any system (0) possesses infinitely many 2-component reductions of the form (|2"j) 
parametrized by two arbitrary functions of a single argument (see examples below). 

On the contrary, the requirement of the existence of nontrivial 3-component reductions 
is already sufficiently restrictive and, in particular, implies the existence of n-component 
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reductions for arbitrary n. This phenomenon is similar to the well-known three-soliton 
condition in the Hirota bilinear approach [231211 (recall that two-soliton solutions exist 
for arbitrary PDEs transformable to Hirota's bilinear form and, therefore, cannot detect 
the integrability) , and the condition of three-dimensional consistency in the classification 
of discrete integrable systems on quad-graphs [I]. One can show, by analyzing equations 
© , © ) that the maximum number of n-component reductions the system © may possess 
is parametrized, modulo changes of variables R l — > f l {R l ) 1 by n arbitrary functions of a 
single argument (notice that this number does not depend on m). Therefore, we propose 
the following 

Definition. A (2 + 1)- dimensional quasilinear system is said to be integrable if it 
possesses n-component reductions of the form ^) parametrized by n arbitrary functions 
of a single argument. 

Remark 2. This definition of integrability implies that the matrices A and B in © 
are not commuting (if we want the system © to be nonlinear and coupled). Therefore, our 
definition is by no means exhaustive: there exist examples of multidimensional quasilinear 
systems which are integrable (in the inverse scattering sense) although do not possess 
n-component reductions. In all these examples the corresponding matrices A and B 
commute: [A, B] = 0. 

Remark 3. The anzatz somewhat similar to the key element of our construction 
was used for finding formal solutions of nonlinear evolution equations in the form of 
exponential series, see (22| and references therein, where solutions were sought as functions 
of real exponentials solving the linearized equation. The possibility to distinguish between 
integrable and nonintegrable equations in the framework of this approach was pointed out 
in |2Zj. 

In section 2 we discuss explicit examples which demonstrate that the above definition 
is indeed very effective for detecting the integrability. 

In section 3 we classify integrable systems of conservation laws in Godunov's form [T5] . 

v t + {fv)x + {g v ) y = 0, w t + {f w ) x + (g w ) y = 0; 
notice that systems of this type automatically possess one extra convex quadratic entropy 

^(v 2 + w 2 ) t + (vf v + wf w - f) x + (vg v + wg w - g) y = 0. 

The integrability conditions constitute a complicated overdetermined system (}4*T]) of fourth 
order PDEs for the potentials / and g. The analysis of this system leads to two possibil- 
ities. 

Quadratic case. There exists a linear combination of / and g which is quadratic in v, w. 
Without any loss of generality one can assume that g is quadratic, say, g = (v 2 — w 2 )/2 
(one has a freedom of Euclidean isometries of the (v, w)-plane to bring g to the canonical 
form). In this case our equations take the form 

v t + (fv)x + v y = 0, w t + (f w ) x -w y = 0; (10) 
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the corresponding integrability conditions reduce to the system (|44|) of fourth order PDEs 
for the potential / which can be solved explicitly (see Sect. 3). Notice that in the new 
independent variables £ = — (t + y)/2, r\ = — (t — y)/2 the system (fTUj) takes the form 



(fv)x = v& {f w ) a 



w 



'/• 



which is manifestly Hamiltonian in the new variables V — f v , W = f w : 

V x = (F v ) ( , W x = (F w ) v . 

Here F is the Legendre transform of /, F = vf v + wf w — f. Thus, we obtain a complete 
description of the class of integrable two-component (2+l)-dimensional Hamiltonian sys- 
tems of hydrodynamic type. An independent treatment of the Hamiltonian case is given 
in Sect. 4 where, in particular, the integrability conditions in terms of the Hamiltonian 
density F are derived (formulae ([56)1 ). 

Harmonic case. Here both / and g are harmonic functions. Further analysis shows that 
there exists a unique system of this type, with / = Re(z In z — z), g = Im(z In z — z), z = 
v + iiv. The corresponding equations are 

VV X + WW X VWy — WVy WV X — VW X VVy + WWy 

vt+ X y~ 2 + f- 2 y =0, w t + x 2 X + y 2 y =0, 

or, in conservative form, 

v t + (In \Jv 2 + w 2 ) x + (arctg-)y = 0, w t — (arctg—) x + (In \Jv 2 + w 2 ) y = 0. 

Remarkably, this two-component system is a disguised form of the nonlinear wave equa- 
tion. To see this we change to polar coordinates v = rcos#, w = rsin#: 

By = — - + r sin6 9 t — cos9 r t , 6 X = — + r cos 6 Q t + sin 6 r t ; (11) 

these equations lead, upon cross-differentiation, to the nonlinear wave equation (r 2 ) tt = 
A ln(r 2 ), known also as the Boyer-Finley equation [3]. Solving equations (jlljl for r x and r y 
and cross-differentiating, one obtains another second order equation (r 2 6 t )t = A# which 
can be viewed as a linear wave equation for 9. It is not clear whether the wave equation 
with a general nonlinearity can be written as a two-component first order system; however, 
as demonstrated in it always possesses a simple three-component representation. 
Hydrodynamic reductions of the Boyer-Finley equation were recently studied in ^3], see 
also Sect. 3. 

systems possessing infinitely many hydrodynamic reductions is provided in section 4. 
The corresponding Hamiltonian densities are solution space is 10-dimensional. 



2 Examples 

In this section we list some of the known examples of (2 + 1) -dimensional systems of 
hydrodynamic type which are integrable in the above sense. All these examples turn out 
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to be conservative and, moreover, possess exactly one 'extra' conservation law which is 
the necessary ingredient of the theory of weak solutions. This makes the systems below 
a possible venue for developing and testing mathematical theory (existence, uniqueness, 
weak solutions, etc.) of multidimensional conservation laws which, currently, remain terra 
incognita [Hj. 

Example 1. The dispersionless KP equation, 



(u t - uu x ) x = u 



plays an important role in nonlinear acoustics and differential geometry. Introducing the 
potential u = (p x we obtain the equation cp xt — (f x (f xx = (p yy which takes the quasilinear 
form 

Vy = W X , Wy =V t ~ VV X (12) 

in the variables v = <p x , w = ip y . Looking for reductions v = v (R 1 , R n ), w = 
wlR 1 , ...,R n ), where the Riemann invariants R l satisfy (j2J), one readily obtains 

d iW = fi'diV, A i = v + {ii 1 ) 2 . (13) 

The compatibility condition didjW = djdiW implies 

didjV = -^-diV + -p^-rdjV, (14) 

while the commutativity condition Q results in 

V = (15) 
The substitution of (fTKj) into (|T4*jl implies the Gibbons-Tsarev system for v(R) and fJ, l (R), 

V = -^-r. a& = 2*!%, (16) 

i 7^ j, which was first derived in fT7\ ITH] in the theory of hydrodynamic reductions of 
Benney's moment equations, see also |2T| l29| 180] for further discussion. For any solution 
fi l ,v of the system (fTT)J) one can reconstruct A* and w by virtue of piljl . In the two- 
component case this system takes the form 

q i d 2 v 2 div d 1 vd 2 v 

d ^ = d ^ = dld2V = W^f {1) 

The general solution of this system is parametrized by four arbitrary functions of a single 
argument, indeed, one can arbitrarily prescribe the Goursat data ^(-R 1 ), ^(R 1 ) and 
v(R 2 ), fi 2 (R 2 ) on the characteristics R 2 = and R 1 = 0, respectively. Moreover, the 
system (1 17 J) is invariant under the reparametrizations R 1 — > f 1 (R 1 ), R 2 — > f 2 (R ) where 
J 1 ,/ 2 are arbitrary functions of their arguments. Since reparametrizations of this type 
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do not effect the corresponding solutions, one concludes that two-component reductions 
are parametrized by two arbitrary functions of a single argument. A remarkable feature 
of the system (JT7j) is its multidimensional compatibility, that is, the compatibility of the 
system (jlfij) which is obtained by 'gluing together' several identical copies of the system 
()17j) for each pair of Riemann invariants R l , W . Indeed, calculating dk{djfi l ) by virtue of 
(US)) one obtains 

t djvdkvj^ + //* - V) 

mM (^-AtW-M*)^*-^)' 
which is manifestly symmetric in j, k. Therefore, the first compatibility condition dk{djH % ) = 
dj(dkfi l ) is satisfied. Similarly, the computation of dk(didjv) results in 

_ djvdjvdkV^) 2 + (j£f + ((j*) 2 - /iV - //// - y?\i k ) 

which is totally symmetric in z, j, fc. Therefore, the second compatibility condition dk{didjv) 
9j(didkv) is also satisfied. The general solution of the system (fTBj) depends, modulo trivial 
symmetries R l —>■ f t (R l ), on n arbitrary functions of a single argument. 

We just mention that the system (fT2|) possesses exactly three conservation laws of 
hydro dynamic type: 

Vy = W X , 



w y = v t - (v 2 /2) x , 



(vw) y = (v 2 /2) t + (w 2 /2 - v 3 /3) x . 
Example 2. The Boyer-Finley equation, 

u xy — (e u )tt, 

(notice that the signature here is different from the one discussed in the introduction, 
which makes the analysis much easier) is descriptive of self-dual Einstein spaces with a 
Killing vector Introducing the potential u = <p t , one obtains the equation tp xy = (e^t 
which takes the form 

v t = w x /w, w t = v y (18) 

in the new variables v = tp x , w — e^*. Looking for reductions in the form v = v^R 1 , R n ), 
w = wlR 1 , R n ), where the Riemann invariants R 1 satisfy (j2j), one readily obtains 

d iW = wXdiV, = w(A*) 2 . (19) 

The compatibility condition didjW = djdiW implies 

d A* d 

Mjv = jrrxA v + xbx 9 ^ ( 2 °) 
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while the commutativity condition Q results in 

= ^—d 3 v. (21) 
A-? - A 4 v ' 

The substitution of (J2T|) into (J20J) implies the system for v(R) and A*(i?), 

(A') 2 A ? A'A'fA' • A') 

^ = i-^-^, 0^ = (A ;_+ )2 I ffrdjV, (22) 

which, in a somewhat different form, was first discussed in For any solution \ % ,v of 
the system one can reconstruct /z\ io by virtue of (flUjl. The system (j22j) is compatible, 
with the general solution depending on n arbitrary functions of a single argument (modulo 
trivial symmetries R l — > f l (R 1 )). 

The system fllSjl possesses three conservation laws of hydrodynamic type: 

W t = Vy, 

v t = {hiw) x , 

(vw) t = w x + {v 2 /2) y . 

Example 3. An interesting integrable modification of the Boyer-Finley equation is 
the PDE 

u xy = (d? - cd 2 x )(e^ 

c = const (see |SH ESj for further examples of this type). Introducing the potential 
g u = 'fix, one obtains the equation (p t t — cip xx = (hi(p x ) y which takes the form 

v t = cw x + Wy/w, w t = v x (23) 

in the new variables v = (fit, w — (p x . Looking for reductions in the form v = v{R 1 , R n ), 
w = wlR 1 , ...,R n ), where the Riemann invariants R 1 satisfy (J2J), one obtains 

d i v = X i d i w, = w{{\ i ) 2 -c). (24) 

The compatibility condition didjV = djdiV implies 

dodjW = — ^ — r-doW + — r-dnw, (25) 
while the commutativity condition Q results in 

d s \* = ( f )2 ~ C ^. (26) 

The substitution of (j2Tlj) into (J25J) implies the following system for w(R) and X l (R), 
ai< ^ l ?-c d 3 w (AQ 2 + (A J ) 2 -2c a 

OjA = -r- r : — , OiOjW = —-. — rr: OiWOjW. 27) 

3 V-A l w' 3 w{\i - A 4 ) 2 3 v ' 
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For any solution X l ,w of the system (|27j) one can reconstruct /i\t> by virtue of ([24)1 . The 
system ()27|) is compatible, with the general solution depending on n arbitrary functions 
of a single argument (modulo symmetries R 1 — ► f l (R 1 )). 

The system (12 3 J) possesses three conservation laws of hydrodynamic type: 

v t = cw x + (lnw) y , 
w t = v x , 

(vw) t = (v 2 /2 + w 2 /(2c)) x + w y . 

Notice that when c > the x-fiux of the third conservation law is a convex function of 
the previous two x-ffuxes (convex entropy). This example can be of particular interest for 
the general theory of multidimensional strictly hyperbolic conservation laws. 

Example 4. We are going to demonstrate that the method of hydrodynamic reduc- 
tions is in fact the effective classification criterion. As an illustration of our approach we 
will classify integrable nonlinear wave equations of the form 

U X y = (/(«))«; 

it will follow that f(u) = e u is the only nontrivial possibility. Introducing the potential 
u = <ft, one obtains the equation ip xy = (f(<ft))t which takes the form 

vt = w x , f'(w)w t = v y (28) 

in the new variables v = <p x , w — <ff Looking for reductions in the form v = v(R l , R n ), 
w = w{R 1 , R n ) where the Riemann invariants R 1 satisfy (J2J) one readily obtains 

d i v = d i w/\\ J = f'(w)(X) 2 . (29) 

The compatibility condition didjV = djdiV implies 

3 {X 3 — X l )X l (X % - X3)X3 3 y ' 
while the commutativity condition Q results in 

f" ( X i ) 2 

= jjnh 9 ^ < 31 > 

The substitution of (j2"Tj) into (|20|l implies the system for w(R) and X l (R), 

f/i (X^ 2 f" X i X-' 
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A direct computation of dk{dj\ l ) implies 

' f"\ 2 (A') 2 (A'((A^) 2 + (A fc ) 2 )+AJA fc (A-'+A fc )-4A'AJA fc ) a o 
/' J (A3-A fc ) 2 (Ai-A i )(A fc -A l ) OjWOkW. 

The compatibility condition d^{dj\ l ) = dj{d^X % ) is equivalent to the requirement that the 
above expression is symmetric in j, k. Since the second term is manifestly symmetric, one 
has to require 

7"Y 
7) ^° 

to ensure the compatibility This implies / = ae bu + c, which is essentially the case of 
Example 2. 

Example 5. There exist remarkable examples with a fairly simple structure of hy- 
drodynamic reductions. Let us consider the system [HI] 

v t =w x , w t = wv x - vw x + Vy. (33) 

Looking for reductions in the form v = v(R l , ...,R n ), w = wlR 1 , ...,R n ) where the Rie- 
mann invariants R l satisfy (J2J) one obtains 

diiv = XdiV, fi* = (A*) 2 + v\ l - w, 

so that 

d A* d 

didjV = r-OiV + — - — 7-rdjV. 

3 \i - X 1 X 1 — \i 3 
The commutativity condition @ implies 

djX 1 = —djV, didp = 0. 

Hence, 

v = j2f k (R% \ l = <p\R l )-J2f k (R k )i 

k k 

where f l {R l ) and (p l (R l ) are arbitrary functions of a single argument. Further properties of 
these reductions were investigated in jM]. We just mention that the system possesses 
three conservation laws of hydrodynamic type: 

v t = w x , 

(w + V 2 ) t = (vw) x + Vy, 

(2vw + v 3 ) t = (v 2 w + w 2 ) x + {v 2 ) y . 
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Example 6. The method of hydrodynamic reductions carries over to multicomponent 
situation in a straightforward way. Here we give details of calculations for the 3-component 
system first proposed in [35], see also fU\ : 



a t + (av) x = 0, v t + vv x + w x = 0, w y + a x = 0. (34) 

Looking for reductions in the form a = a(R l , R n ), b = b(R l , R n ), w = w(R l , R n ) 
where the Riemann invariants R l satisfy (J2J one obtains the relations 

diW = -(X i + v)diV, did = /2 i (X i + v)diV, /J = - ° . (35) 
The compatibility condition didjW = djdiW implies 

d&v = jj^div + x^fyv, (36) 
while the commutativity condition (JHJ) results in 

The substitution of (|S7|) into (|3l)|) implies the system for v(R) and X l (R), 

One can verify that the remaining compatibility conditions didja = djdid are satisfied 
identically. For any solution A*, v of the system ()38j1 one can reconstruct w, a, fi 1 by virtue 
of ()35j) . The system (|3*H|) is compatible, with the general solution depending on n arbitrary 
functions of a single argument (modulo symmetries R l — > f l (R 1 )). 

Notice that the system (jMj) possesses four conservation laws of hydrodynamic type: 



a t + (av) x = 0, 
v t + (v 2 /2 + w) x = 0, 
w y + a x = 0, 

(aw + av 2 ) x + (w 2 /2) y + (au) t = 0. 

Equations (}3~l"j) can be generalized as follows: 

a t + (av) x = 0, u t + + w x = 0, w y + p(a) x = 0. 

(so that one obtains isentropic gas dynamics in the limit x = y). It can be shown that 
this system passes the integrability test if and only if p" = 0, which leads to ()34jl . 
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3 Classification of integrable systems of conservation 
laws with a convex quadratic entropy 

In this section we discuss systems of conservation laws in Godunov's form 

v t + (fv)x + {g v )y = 0, w t + (f w ) x + {g w ) y = 0; (39) 

here f(v,w) and g(v,w) are given potentials. Systems of this type automatically possess 
one extra convex quadratic entropy 

i(> 2 + w 2 ) t + (vf v + wf w - f) x + (vg v + wg w - g) y = 0. 

Equations ()39|) can be written in the matrix form (^) with 



u ! V ) , A = ( fvv fvw ), B = ( 9vv 9 

W I \ J vw J ww J \ 9vw 9 



vw 

WW 



in what follows we assume that the commutator 

[A, B] = 



s 
-s 

is nonzero, that is, s = g vw (fvv — fww) — fvwidw — 9ww) 7^ 0; otherwise, the system possesses 
no nontrivial n-component hydrodynamic reductions. The integrability conditions lead 
to an overdetermined system of fourth order PDEs for / and g. A careful analysis of this 
system shows that there exist two essentially different possibilities: 

- g is quadratic in v , w (quadratic case) or 

- both / and g are harmonic functions (harmonic case). 

Remarkably, in both cases the equations for / and g can be solved in a closed form. 
We hope that these examples would provide a good venue for developing and testing 
the general theory of multidimensional conservation laws (breakdown of solutions, weak 
solutions, etc). 

The integrability conditions can be derived in the standard way. Looking for reductions 
of the system (j2*3j) in the form v = v(R 1 , ...,R n ), w = w{R}, ...,R n ) where the Riemann 
invariants satisfy equations (J2J), and substituting into (|23|) . one arrives at 

{X + fvv+^gw)diV + (f V w+^gvw)diW = 0, (fvw+V t gvw)diV + (\ t + fww+^gww)d i w = 0, 

(no summation!) so that A* and // satisfy the dispersion relation 

(A* + fvv + A^toXA* + fww + gww) (fvw fj> gvw) ■ 
Setting d{V = ip l diW one obtains the following expressions for A* and \x l in terms of (p l , 

(fvvQvw fvwgvv)i}P ) (fvvgww fwwgvv)^P (fvwgww fwwgvw) 



A' 



g vw {\ - (y?*) 2 ) + (g vv - gwwW 

f vw (l - (y') 2 ) + (fvv - fww)^ 

SW>(1 - (v 9 *) 2 ) + (9w - 9wwW 
12 



which define a rational parametrization of the dispersion relation, 
conditions of the equations d{V = f l diW imply 



didjiv 



— :OiW H : rOiW 



If J — f l 



f — If J 



The compatibility 



(40) 



while the commutativity conditions Q lead to the expressions for djip 1 in the form djip 1 = 
(...)djW. Here dots denote a rational expression in if 1 , if* whose coefficients are functions of 
the second and third derivatives of / and g. We do not write them down explicitly due to 
their complexity. To manipulate with these expressions we used symbolic computations. 
One can see that the compatibility condition dkdjf l —djdkf % = is of the form PdjwdkW = 
0, where P is a complicated rational expression in ip l , f?, f k whose coefficients depend on 
partial derivatives of / and g up to fourth order. Requiring that P vanishes identically we 
obtain the overdetermined system of fourth order PDEs for / and g. This system yields 
the following expressions for the fourth derivatives of /: 

sfvvvv fvvv\P'\9ww 9vv) fvvw ^xfvv fww)gvvw ~\~ 

^dvwi^fvvv ^fvww) ^fvw9vww\ 

9vvv\\f WW fvv)f vvw 2 fvw ( fvww fvvv)\ ~\~ 

6 fvvw(fvw9vvw 9vw fvvw) 1 
sfvvvw fvvw\fvw9vvv ~t" < ^ > 9vw fvww ^fvw9vww) fvvv\P'\fww fvv)9vww 

2 \9ww 9w) fvww ~\~ fvw{^9www ~\~ 9vvw) 9vw(^fwww ^fvvw )] > 
sfvvww (fi'uiui 9vv) (fvww fvvw fvvvfwww) ~\~ 

\fvv fww)(fvww9vvw 9vvvfwww) (^"0 

^9vw{fvvw fvww) ^ fvw {fvww9vww fvvw gvvw) ■> 
sfvwww fvww(fvw9www 3g vw fvvw <$f vw g vvw ) ~\~ fwww\^\fvv fww)9vvw 

2(5Wi 9ww) fvvw fvw\9vvv 9vww) 9vw(fvvv 2 fvww)\i 
sfwwww fwww [2 i^9vv gww) fvww ~\~ ^\fww fvv)gvww ~\~ 

^9vw(fwww ^fvvw) 2 fvw9vvw\ 

9www\f\fvv fww) fvww 2 fvw ( fvvw fwww)\ 

Q fvww(^fvw9vww gvw fvww ) ) 

here s = g vu) {f vv — f ww ) — fvw(gw — 9ww) ^ 0. Notice that equations (j4TJ) 4 and (JUJs can 
be obtained from (|41jlo and (|41jW by interchanging v and w. Analogous expressions for 
the fourth derivatives of g can be obtained from (|41|) by interchanging / and g. Moreover, 
one has five quadratic relations among the third derivatives of / and g: 



fvwg 
fvw g 

fvw 9 



www 
www 
www 



f 



u !uw9vvv fvww9www 

fwww gvw fvw gvvw 

f www9vvv fvvw 9vww 

fvvv9vww fvww9vvv ~ 

fvvw9www fwww9vvw 



9v 



I wwwyvww 



f 
f 

fvww9vvw 

o, 

= 0. 



= 0, 

r ow gvw o, 

0, 



(42) 
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This system of PDEs for / and g is not in involution; a careful analysis of the quadratic 
relations (|42|) leads to two essentially different possibilities: 

Quadratic case. The third derivatives of g are proportional to the corresponding third 
derivatives of /, 

Qvvv fvvv 5 Qvvw f^f vvw) gvww f^fvwwj gwww fJ'fwww 

Substituting these relations into the remaining integrability conditions one can show that 
\i must be constant. Therefore, g — /// is at most quadratic in v,w. Without any loss of 
generality one can assume that, say, g is quadratic. 
Harmonic case. Here 

fvvv fvww 0, fwww fvvw 0, gvvv gvww 0, gwww gvvw 0, 

which imply that Af and Ag are constants. Without any loss of generality one can 
assume that both / and g are harmonic, Af = Ag = 0. 

Remark 4. There exists an obvious group of equivalence transformations which 
preserve the integrability and leave equations (ffilj) form- invariant. These are, first of all, 
orthogonal transformations of the (v, ty)-plane, generated by translations and rotations. 
Secondly, these are linear changes of the independent variables in (jHSj) , 

x — > aux + a 12 y + a^t, y — > a 2 \x + a 22 y + a 23 t, 

which induce the transformations 

v 2 + w 2 , v 2 + w 2 

J -> ant + a i29 + «i3 ^ , 9 -> a 21 f + a 22 g + a 23 • 



The classification below is carried out up to this natural equivalence. 
3.1 Quadratic case 

Setting g(v,w) = av 2 + 2bvw + cw 2 and introducing C = b(f ww — f vv ) + (a — c)f vw one 
can rewrite equations (|41jl as follows: 

C fvvw ^^{^fvvw ^ fvvv fvww fvvv) 2(o c) fvvv fvvw j 

C fvww b(3 fvvw fvww 2 fvvv fvvw fvvv fwww ) "I - 2(fl c) fwv fvwwi 

C fwww 26(/ tIUIlu fvvw) C)(fvvwfvww fvwfwww)i (^3) 

C fvwww b(J$ fvww fvvw 2 fwww fvww fvvv fwww) 2{o c) fvvw fwww j 

C fvjwww ^O^fvww 2 fwww fvvw fwww) 2(o c) fvww fwww ; 

(the corresponding equations for g are satisfied identically). For any a, 6, c this system 
is in involution and its solution space is 10- dimensional. Indeed, the values of /, its 
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first, second and third derivatives can be choosen arbitrarily, while the fourth and higher 
derivatives are determined by virtue of (|43|). Diagonalizing the quadratic form g by a 
linear orthogonal change of variables v,w, one can set 6 = 0. In this case equations (|43J1 
simplify to 

fvwf vvvv ^fvwfwwi 
fvwfvww ^fvwfvwwi 
fvwfvvww fvvwfvww fvvvfwwwi (^4) 
fvwfvwww 2 f vvw fyj ww , 
fvwfwwww ^fvwwfwww 

The first two equations imply that fvw/ fvw = cons t- Similarly, the last two equations 
imply fwwwl ftw = const. Setting f vw = e one can parametrise the third derivatives of / 
as follows: ^ ^ 

fvw —TT16 , fww fvww fwww ' (4^) 

here m, n are arbitrary constants. The compatibility conditions of these equations plus 
the equation ()44jh result in the following overdetermined system for e: 

It is worth mentioning that the system ()46|) arises in a completely different context in 
projective differential geometry constituting the projective Gauss-Codazzi equations of 
the Roman surface of Steiner which is known as the only quartic in P 3 containing a two- 
parameter family of conies, see e.g. ^2]- At the moment we have no explanation of this 
remarkable coincidence. 

Solving the first (Liouville) equation for e in the form 



2 



4 p'{y)q'{w) 



ran (p(v) + q(w)) 2 

and setting 

(pf/ 2 = ^iP(p), (q') 3/2 = V^Q(q), (47) 

(here P(p) and Q{q) are functions to be determined), one obtains from the last two 
equations (}4T)j) the following functional-differential equations for P and Q: 

P"(p+q) 2 -4P'(p+q)+6P = 2Q'(p+q)-QQ, Q"(p+q) 2 -4Q'(p+q)+QQ = 2P'(p+q)-QP; 

these equations imply that both P and Q are cubic polynomials in p and q, 

P = ap 3 + bp 2 + cp + d, Q = aq 3 — bq 2 + cq — d, 
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where a, b, c, d are arbitrary constants. In the general case, using translations and scalings, 
equations (}4Tj) can be brought to the form (y') 3 = (y 3 — 3Xy 2 + 3y) 2 and solved in terms 
of elliptic functions |2j: 



V = ~x o ,/ » \ ' #3 



4-3A 2 



X-3p'(z;0,g 3 y 27 
In the simplest case m = n = equations (|46|) imply 

e = (av + P)( r yw + 5), 
and the elementary integration of (|45jl results in 

\ «7 2 2 «^ 2 2 nc 

j [v, w) = —v w + —v w + —vw + povw; 



here a, (3,j,5 are arbitrary constants. Using the equivalence transformations one can 
reduce / to either / = v 2 w 2 (if both a and 7 are nonzero) or / = vw 2 (if a = 0). The 
corresponding equations take the form 

u t + 2(W% + w y = 0, iu t + 2(u 2 u>) x - tOj = 

and 

v t + {w 2 ) x + v y = 0, iu t + 2(vw) x -w y = 0, 

respectively (in both cases g = (v 2 — w 2 )/2). 
If m = 0, n^O, equations (J46|) imply 



e = (cru + P)<p'(w) 

where <p"' = an((p') 2 , a, (3 — const. Therefore, (p(w) = — -^((w;0,g 3 ) + 7, where C( z ) 
—p(z) and #3,7 are constants. The elementary integration of the equations (|4"5j) gives 

f = ^(av + P) 2 v(w). 
This reduces to the previous case if n = 0. 

3.2 Harmonic case 

We have f vv = -f ww , g vv = -g ww , so that 

fvvv fvwwt fwww fvvwi gvvv gvwwt gwww g 



WW ■ 



Differentiating these relations by virtue of (|4~T|) (and the analogous equations for g) one 
arrives at the additional constraints 

fvvPl + fvwP2 + g ww P3 = 0, - f vw Pl + f vv p 2 + g V wP3 = 0, 

-g V wPi + g ww p2 + fvwPA = 0, g ww p\ + g vu ,P2 + f m p± = 
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where 

Pi fvww 9 vww fww 9 WW J P2 fvww9vvw fvvw9vwwi 

fww fvww ' P^ 9vvw~^~9 



Solving the linear system for pi (notice that the corresponding 4x4 matrix has rank 
three) one obtains 



Pi fvww9vww fvvw9vvw (fvw9vw fvv9ww) 1 

P2 fvww 9 vvw fvvw9vww y> (fvv9vw fvw9ww, 

D , = f 2 + f 2 = u 2 (f 2 + f 2 ) 
fo J vvw 1 J vww \J vw 1 J vv J ' 

9 vvw 9 vww \9vw 9ww)' 

Introducing the two-component vectors 

. fvww \ ( 9vww \ ( fvw 

\ „ 9 

e l — \ J- I j e 2 — „ , s l — x , S 2 



(4* 



vw 

WW 



j toui / \ yww / \ J ww j \ y-u 

one can rewrite (|48|) in vector notation as follows 

(ei,ei) = /i 2 (si,si), (e 2 ,e 2 ) = /i 2 (s 2 ,s 2 ), (e x , e 2 ) = //(si, s 2 ), e x A e 2 = -/i 2 Si A s 2 

where (, ) is the standard Euclidean scalar product. Therefore, the vectors e±, e 2 and si, s 2 
are related by a composition of a scaling, rotation and reflection, that is, 

s 2 





' X 




' X 




1=1 




_x ) Sl ' 62 = 1 







or, explicitly, 

fvWW fvw Y fyjWt fwW ^ f VW f'l 



(49) 

9 vww 9vw Y 9wwi 9 ww Y 9vw X 9v 



WW 1 
WW 



(notice that equations for / and g coincide). The compatibility conditions of these equa- 
tions imply the following equations for X and Y, 



X v — 2XY, X w — X — Y , Y v — Y — X , Y w — 2.XY. 
whose general solution (up to translations in v and w) is 

X = - „ W Y 



V 2 + W 2 ' V 2 + w 2 

Substituting these expressions into (|4T)|) and integrating the corresponding linear system 
for the harmonic function / (this integration simplifies if one changes to the complex 
variables z = v + iw, z — v — iw) one readily obtains that / must be a linear combination 
of the real and imaginary parts of the function z In z — z. The same result holds for g. 
Therefore, without any loss of generality one can set / = Re(z In z— z), g = Im(z In z— z). 
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3.3 Reductions of the system ( II lj) 

Here we briefly discuss reductions of the system , 

9 y — — - + r sin 9 9 t — cos 9 r t , 9 X = — + r cos 6 9 t + sin 9 r t . 

Assuming 9 = 9(R 1 , R n ), r = r(R l , R n ) where the Riemann invariants satisfy the 
equations 

R l = A l (-R) R v Rl y = ^{R) R\ 

and substituting into (fTTj) one arrives at 

(p* - r sin 9)dt9 + (\'/r + cos 9) dp = 0, (A' - r cos 9)8$ - (tf/r + sin 0)d*r = 0, 

(no summation!) Hence, the characteristic speeds \ l and \i % satisfy the dispersion relation 
which takes a particularly simple form 

(A*) 2 + (//) 2 = r 2 . 

Parametrising A* and fi l in the form X 1 = r cos if 1 , fi l = r sirup 1 we obtain 

cos 9 + cos tp i 

di9 = —q : T diU ' 50 

sin 9 — sm if 1 

U = lnr. The compatibility conditions of these equations together with the commutativ- 
ity conditions (JHJ imply the system 

WU = --rT$L;, d j <p* = cat ( £^-d j U (51) 
sin v * L 

which is a trigonometric version of the Gibbons-Tsarev system ()16|) . As shown in ^3] 
this system is, in fact, equivalent to (|16p. and its solutions can be constructed from the 
known solutions of the Gibbons-Tsarev system ^3 EE]- Once the solution of (|51j) is 
known, the corresponding polar angle ^ can be calculated from the equations (|5()|) which 
are compatible by construction. 



4 Classification of integrable Hamiltonian systems of 
hydrodynamic type in 2 + 1 dimensions 

In this section we classify Hamiltonian systems 

v t = (K)x, w t = (h w ) y (52) 

which possess infinitely many hydrodynamic reductions. Here h(v,w) is the Hamiltonian 
density. Notice that any system of the form (}52|) possesses one extra conservation law 

ht = (hl/2) X + (hl/2)y. 
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As mentioned in the introduction, Hamiltonian systems are related to the quadratic case 
of Sect. 3 by virtue of the Legendre transform. However, we find it instructive to treat 
the Hamiltonian case independently to better illustrate our method. 

Looking for reductions in the form v = v (R 1 , R n ), w = u^i? 1 , R n ) where the 
Riemann invariants satisfy the equations 

R! X = X\R)R\ ) Rl = ^(R)Rl (53) 

and substituting into (JS2J) one arrives at the equations 

(1 - \ l h vv )diV = X l h vw diW, (1 - fi l h ww )diW = fi l h vw diV, 

(no summation!) so that A 1 and \i l satisfy the dispersion relation 

(1 - Xh vv ){l - ^h ww ) = Ay^L- 

We require that the dispersion relation is nondegenerate (as a conic), that is, h vw ^ 
0, h% w — h vv h ww 7^ 0. Setting d(v = (p l diW, we obtain the following expressions for A 1 and 
\x % in terms of (p l 



LD l ■ 1 

A 1 = i F-n~ . A** 



"vw 



which define a rational parametrization of the dispersion relation. The compatibility 
conditions of the equations d\v = (p l diW imply 

didjW = -^f- — T diW H — — :djW (54) 

while the commutativity equations (jSJ) lead to the following complicated expressions for 

P) , n i _ ^(hvw+V 3 h vv )(h ww +ip z h vm )djW r/, , ir \ j . i _!./„»/, 1-1- 

U 3 { T ~ h vw (h% w -h vv fH uw )(<p i -<pi) Y\""vvw T <+> 'hivvjf T »«» T ^ "-wn«J T 

(55) 

(h vw +ip l h vv )(h ww +ip3h vw )djW r/, j, \ »• , j, -I 

h vw {hl w -h vv h ww )(^—pr) lV L vww T V (Ivvw)V T Uwww T <f ^totoJ • 

One can see that the compatibility condition dkdjip l —djdk<fi l = is of the form PdjwdkW = 
where P is a rational expression in (f 1 , ip 3 , (p k whose coefficients depend on partial deriva- 
tives of the Hamiltonian density h(v,w) up to fourth order. Requiring that P vanishes 
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we obtain the overdetermined system of fourth order PDEs for the density h: 

7, 7,2 _ 97, 7, /, /, _ 7, 7, 7,2 

\'-JI"VV' l 'VW' b vvW ^'l"VV' <"WW' "VVV 1 b vvw ' b vw' b ww' b vwi 



h (h 2 -h h )h = -h h (h h + h h ) 

' "vwy "vw "vv ' "ww J ' "vww ' "vw ' "vvv \ 1 "vv '"www 1 ' "ww ' "vvw ) 

j-lh 2 7)2 — 97) 7) 7) 7) -I- 7)2 7) h 

~ vw vvw vv ' "ww' "vvv ' "vww ~ '"vw VVV "'VWW} 

h (h 2 -h h )h = Ah 2 h h 

'"vw y^vw vv ' "ww ) ' "www ^"vw h vvw' L vww 

hvvhvvwy^vw^www ~t~ ^ww^vww) h , ww^ , wv( y ^ , vw^ , vww ~t~ h , wh , www^)i 

h (h 2 -h h )h — —h h (h h +h h ) 

' "vw \ ' "vw ' "vv' "ww J ' "vwww '"vw'"www\'"ww'"vvv 1 ' "vv ' "vww J 

-Lq/,2 /,2 _ 97, 7, 7, 7, 1 7,2 u u 

~ vw vww vv 1 "ww' "www' "vvw ~ ' "vw viww' "VVW ) 



h (h 2 -h h )h = 47) 7) (h h - h h ) 

IL VW\' L VW ' L VV' L WW )' L WWWW ^' "VW' "WWW\' "VW' L vww ' L ww' L vvw) 

j-ST) 7) 7)2 — 97) 7) 7) h — h h h 2 

\ u '"ww'"vw'"vww "vv' "ww' "www' "vww '"vw'"w'"www 



(56) 



It was verified that this system is in involution and its solution space is 10-dimensional. 
As mentioned in the introduction, the Legendre transform identifies the systems (|56)l and 
(PI). 



5 Concluding remarks 

We have demonstrated that the existence of hydrodynamic reductions describing nonlinear 
interactions of n > 3 planar simple waves can be viewed as the effective integrability 
criterion. The most natural problems arising in this context are the following: 

1. Classify multicomponent (2 + l)-dimensional integrable quasilinear systems. The main 
difference from the two-component case is that the dispersion relation (|7J) will no longer 
define a rational curve. 

2. The recent publication ^3] suggests that the method of hydrodynamic reductions 
carries over to 3 + 1 dimensions. It would be extremely interesting to obtain further 
examples (classification results) of dispersionless integrable systems in many dimensions. 

3. Nonlinear interactions of n simple waves can be viewed as a natural dispersionless 
analogue of 'n-gap' solutions. It would be desirable to obtain an alternative description 
of these solutions as the 'stationary points' of the appropriate 'higher symmetries'. 

We hope to address these questions elsewhere. 



Acknowledgements 

EVF is grateful to A. Fokas for drawing his attention to the problem of integrability 
of multi-dimensional quasilinear systems. We also thank C. Dafermos, M. Dunajski, J. 
Gibbons, B. Keyfitz, M. Pavlov, D. Serre and C. Trivisa for their interest and clarifying 
discussions. This research was initiated when both authors were attending the research 



20 



semester on Hyperbolic Conservation Laws at the Isaak Newton Institute in the spring of 
2003. It is a great pleasure to thank P. LeFloch for the invitation to participate in this 
program. 



References 

V.E. Adler, A.I. Bobenko and Yu.B. Suris, Classification of integrable equations on 
quad-graphs. The consistency approach, |nlin .SI/0202024 

M. Abramowitz and LA. Stegun, Handbook of mathematical functions, Nauka, 
Moscow, 1979. 

M. Blaszak, B.M. Szablikowski, Classical R-matrix theory of dispersionless systems: 



II. (2+l)-dimension theory, nlin.SI/0211018 



CP. Boyer and J.D. Finley, III Killing vectors in self-dual, Euclidean Einstein spaces. 
J. Math. Phys. 23 (1982) 1126-1130. 

M. Burnat, The method of Riemann invariants for multi-dimensional nonelliptic sys- 
tem. Bull. Acad. Polon. Sci. Sr. Sci. Tech. 17 (1969) 1019-1026. 

M. Burnat, The method of Riemann invariants and its applications to the theory of 
plasticity. I, II. Arch. Mech. (Arch. Mech. Stos.) 23 (1971), 817-838; ibid. 24 (1972), 
3-26. 

M. Burnat, The method of characteristics and Riemann's invariants for multidimen- 
sional hyperbolic systems. (Russian) Sibirsk. Mat. Z. 11 (1970) 279-309. 

C. Dafermos, Hyperbolic conservation laws in continuum physics, Springer- Verlag, 
2000. 

L. Dinu, Some remarks concerning the Riemann invariance, Burnat-Peradzyhski and 
Martin approaches, Rev. Roumaine Math. Pures Appl. 35, N3 (1990) 203-234. 

B.A. Dubrovin and S.P. Novikov, Hydrodynamics of weakly deformed soliton lattices: 
differential geometry and Hamiltonian theory, Russian Math. Surveys, 44 , N6 (1989) 
35-124. 

B.A. Dubrovin, Geometry of 2D topological field theories, Lect. Notes in Math. 1620, 
Springer- Verlag (1996) 120-348. 

E.V. Ferapontov, Integrable systems in projective differential geometry, Kyushu J. 
Math. 54, Nl (2000) 183-215. 

E.V. Ferapontov, D. A. Korotkin and V.A. Shramchenko, Boyer-Finley equation and 
systems of hydrodynamic type, Class. Quantum Grav. 19 (2002) L1-L6. 



21 



[14] E.V. Ferapontov and M.V. Pavlov, Hydrodynamic reductions of the heavenly equa- 
tion, Class. Quantum Grav. 20 (2003) 1-13. 

[15] J. Gibbons, Collisionless Boltzmann equations and integrable moment equations, 
Physica 3D (1981) 503-511. 

[16] J. Gibbons and Y. Kodama, A method for solving the dispersionless KP hierarchy 
and its exact solutions. II. Phys. Lett. A 135, N3 (1989) 167-170. 

[17] J. Gibbons and S. P. Tsarev, Reductions of the Benney equations, Phys. Lett. A 211 
(1996) 19-24. 

[18] J. Gibbons and S. P. Tsarev, Conformal maps and reductions of the Benney equa- 
tions, Phys. Lett. A 258 (1999) 263-271. 

[19] S.K. Godunov, An interesting class of quasi-linear systems, Dokl. Akad. Nauk SSSR 
139 (1961) 521-523. 

[20] A. Grundland and R. Zelazny, Simple waves in quasilinear hyperbolic systems. I, II. 
Riemann invariants for the problem of simple wave interactions. J. Math. Phys. 24, 
N9 (1983) 2305-2328. 

[21] F. Guil, M. Manas and L. Martinez Alonso, On the Whitham Hierarchies: Reductions 
and Hodograph Solutions, |nlin.SI/0209"05T| 

[22] W. Hereman, P.P. Banerjee, A. Korpel, G. Assanto, A. Van Immerzule and A. Meer- 
poel, Exact solitary wave solutions of non-linear evolution and wave equations using 
a direct algebraic method, J. Phys. A: Math. Gen. 19 (1986) 607-628. 

[23] J. Hietarinta, Introduction to the Hirota bilinear method. Integrability of nonlinear 
systems (Pondicherry, 1996), 95-103, Lecture Notes in Phys., 495, Springer, Berlin, 
1997. 

[24] J. Hietarinta, Equations that pass Hirota's three-soliton condition and other tests 
of integrability. Nonlinear evolution equations and dynamical systems (Kolymbari, 
1989), 46-50, Res. Rep. Phys., Springer, Berlin, 1990. 

[25] J. Hietarinta, A search for nonlinear equations passing Hirota's three-soliton condi- 
tion, J. Math. Phys. 28 (1987) 1732-1742, 2094-2101, 2586-2592; 29 (1988) 628-635. 

[26] S. Canic, B.L. Keyfitz and E.H. Kim, Mixed Hyperbolic-Elliptic Systems in Self- 
Similar Flows, Boletim da Sociedade Brasileira de Matematica 32 (2002) 1-23. 

[27] K.R. Khusnutdinova, Exact solutions describing interaction of solitary waves of non- 
integrable equations, in: Nonlinearity and Geometry, PWN, Warsaw (1998) 319-333. 

[28] A. Majda, Compressible fluid flows and systems of conservation laws in several space 
variables, Appl. Math. Sci., Springer- Verlag, NY, 53 (1984). 



22 



M. Manas, L. Martinez Alonso and E. Medina, Reductions and hodograph solutions 
of the dispersionless KP hierarchy, J. Phys. A: Math. Gen. 35 (2002) 401-417. 

Lei Yu, Reductions of dispersionless integrable hierarchies, PhD Thesis, Imperial 
College, London, 2001. 

M. Manas and L. Martinez Alonso, A hodograph transformation which applies to 
the heavenly equation, nlin. SI/0209050 

A.V. Mikhailov and R.I. Yamilov, Towards classification of (2 + l)-dimensional inte- 
grable equations. Integrability conditions. I. J. Phys. A 31, N31 (1998) 6707-6715. 

A.V. Mikhailov, A.B. Shabat, and V.V. Sokolov, The symmetry approach to classi- 
fication of integrable equations, in: What is integrability? Springer Ser. Nonlinear 
Dynam., Springer, Berlin, (1991) 115-184. 

M.V. Pavlov, Integrable hydrodynamic chains, nlin. SI/0301010 

M. V. Pavlov, Classification of the integrable Egorov hydrodynamic chains, submitted 
to Theor. and Math. Phys. (2003). 

Z. Peradzyhski, Riemann invariants for the nonplanar /c-waves. Bull. Acad. Polon. 
Sci. Sr. Sci. Tech. 19 (1971) 717-724. 

Z. Peradzyhski, Nonlinear plane fc-waves and Riemann invariants. Bull. Acad. Polon. 
Sci. Sr. Sci. Tech. 19 (1971) 625-632. 

S.P. Tsarev, Geometry of hamiltonian systems of hydrodynamic type. Generalized 
hodograph method. Izvestija AN USSR Math. 54, N5 (1990) 1048-1068. 

E.V. Zakharov, Dispersionless limit of integrable systems in 2 + 1 dimensions, in 
Singular Limits of Dispersive Waves, Ed. N.M. Ercolani et al., Plenum Press, NY, 
(1994) 165-174. 



23 



